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1. Introduction 

In braneworld scenarios our Universe is pictured as a hyper surf ace living in a higher dimen- 
sional spacetime, with our standard model particles bound to the hypersurface. In string 
theory the natural choice for such hypersurfaces are the D-branes Q, or the boundaries 
of spacetime Q. However, the braneworld idea was first proposed in the context of field 
theory models, using topological defects to construct bound states of scalar fields on their 
world- volume |Q. 

In condensed matter systems, it is well known that fermions can be trapped by impu- 
rities in the system, in the sense that the lowest energy states are localised around these 
impurities (see, for instance Q). In a similar way, topological defects may lead to localised 
fermion states ||5|, |6| and one can construct the matter sector of the standard of the standard 
model on such defects [0]. ^. 

Multiple branes may exist along the extra dimensions, and so it is likely that brane- 
brane collisions will take place. This has prompted various scenarios for the creation of 
the cosmic microwave background anisotropies reheating |]l^] and baryogenesis. In 
the latter case, the baryon asymmetry is created as the localised fermions carried by the 
colliding branes are transferred in a CP-asymmetric way. 

In [|l^], we extended the work of studying collisions of branes represented by topo- 
logical defects in 1-1-1 dimensional (j)'^ theory (a kink), on which fermions were localised. 

^Note, that there are other dynamical ways of ensuring locahsation to the brane, for instance tel. 
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We found that fermions are indeed transferred in kink-anti-kink collisions, but also when 
bouncing a kink off a boundary. These boundaries can themselves be thought of as branes 
||2|. The final state resulting from kink-boundary and kink-kink collisions was rather sen- 
sitive to the collision velocity and the scalar-fermion coupling, with an oscillatory rather 
than monotonic dependence for the number of transferred fermions [12, pd| ]. 

In 12 1, only the zero momentum fermion modes were considered, and we now 
extend the treatment further by including the bound-states of fermions with non-zero 
momenta along the brane dimensions. The system we simulate is that of a domain wall in 
five dimensions colliding with a (parallel) spatial boundary, we are then able to evaluate 
the particle spectrum of fermions resulting from such a collision, as viewed by braneworld 
observers on either the brane or boundary. 



2. Scalar-fermion model in 4+1 dimensions 

We will model the brane collisions in terms of kinks of a real scalar field 4> in 4+1 dimensions. 
Coupled to this field is a Dirac fermion -0 which, as we shall see, allows for a bound state 
localised on the kink; these correspond to the braneworld fields. The action is 



Sbuik = — j dtdzd X 
to which we add a boundary term 

'^'boundary I dt d X 



^d^cPd^'cP - #7^9^ + ^ 1)' - 



(2.1) 



Jdtd^x 



J 2 = 



For a discussion of this term see |11], |15]. We will use the conventions 
rj^" = diag{-l, 1,1, 1,1), {-Ji^-fu} = ^t],,^, ijj = i/j^q, 



(2.2) 



(2.3) 



and we have singled out the x^ = z direction which will be into the bulk, orthogonal to 
the 3-dimensional branes, all of which will be assumed to be parallel to each other. In 4+1 
dimensional spacetime, the fermions are 4-component complex fields, and we employ the 
chiral representation, for which 
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7J = 


, —ia^ 


) 



-«7o7i7273 




(2.4) 



with the Pauli matrices. We can then split the spinors as 
with 




*iJ = ^(1+7^)*, 



= ^(1-72)^- 



(2.5) 



(2.6) 
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To match the conventions used in [^], we write 

V'L = V'2, ipR = itpi. (2.7) 

where ■01,2 each have two complex component^. 

We assume that the scalar field does not depend on the directions x = Xj, j = 1,2,3 
along the brane but only on the orthogonal direction z. In that case, the scalar equation 
of motion reduces to 

4>{t, z) = dlct^it, z)-\ {4>Ht, z) - 1) z), (2.8) 

whereas for the fermions we have 

i'iit,x,z) + a^dj'ipi{t,x,z) + dz'ip2{t,x,z) - g(j){t,z)-ip2it,x,z) = 0, (2.9) 
'tp2{t,x,z)-a^^j^l^2{t,x,z)+^z'^l^l{t,x,z)+g(|){t,z)^pl{t,x,z) = 0. (2.10) 

We have neglected the fermion back-reaction ig{^^) on the scalar field for reasons of 
simplification. Our main interest is in the behaviour of the fermions, and in the setup here 
back-reaction would only lead to a small correction to the kink profile, indeed, the zero 
modes generate no back-reaction. Of course, in the presence of many fermions, such as in 
a high-temperature thermal state, such effects cannot be discarded (see [13| on how to deal 
with this). 

2.1 The kink and fermion bound states 

The static scalar field equation has a kink solution, 

cPk{z) =tanh(^^^y D=^, (2.11) 

here centered at z = zq with width D. We will operate on the negative half-axis z < 0, 
with the kink incident on a boundary at z = 0. 

To solve the fermion equation of motion in the background of the static kink |14, 12|, 
we will use a separable ansatz for the fermion field, and write 

ipiit,x,z) = J^K{t,z)xiit,x), ip2{t,x,z) = 0. (2.12) 

We find that 

dzJ'Kit, z) + g^{t, z)TK{t, z) = 0, (2.13) 

Xi{t,x) + a^djXi{t,x) = Q. (2.14) 

The J-K is a profile function in the z-direction, and is the term responsible for localizing 
the fermion on the domain wall, 



^ l T[gD + \] 1 _ Nk , . 

\ 2D^T[gD][cos\i{z/D)]9D [cosh(z/Z))]5^ ' ^ ' 



^In the simpler case of hM, these were single-component and real. 
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The equation for xi is that of a right-handed chiral fermion. Depending on the value of 
g whole towers of such localised modes exist but we found that even the second mode 
contributes at most ~ 10 percent |11]. We will here restrict our attention to the lowest 
mode ( p.l5 ) only. 

We follow the standard quantization procedure and expand the component xi terms 
of plane waves along the domain wall, 



where we require that 



d^p 1 
(27r)3 2^p r'- 



bip)e'P''U{p) + d\p)e~'P''V{p) 



p = ujp > \ujp 



\P\- 



(2.16) 



(2.17) 



The first and second terms in the expansion correspond to positive and negative energy 
modes respectively. Substituting this into the equation of motion for xi we find, 



b{-p)[l +p.a]U{-p)e 



—iujt 



d\p)[l-p.q\V{p)t 



iuit 



0, 



with p = p/ujp. The two terms must vanish independently, and so 

p.aVip) = Vip), p.aU{-p) = -U{-p), 

from which we see that 

Vip) = U{p). 
So for instance if we have p = (0, 0, 1), then 



U{p) oc 

We have therefore that 



/ > 0, 



d^p 1 



U{p) oc 







p° < 0. 



6(p)e*P^C/(p) + d\p)e-'P''U{p) 



(27r)3 2ujp 
where we choose the normalization 

U\p)U{p) = lujp. 

This last relation fixes the normalization of the operators h and d so that, given 

^{t,x,z) = TKit.z) I . ° I , / dzT*K(t,z)TK{t,z) = 1, 

the canonical commutation relations 

{*„(t,x,z),^'j3(t,^',/)} = bo,f3b{x-^)b{z- z'), 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



(2.25) 
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reduce to 

{Xia{t,x),x\p{t,x')} = Sap5{x-x!), (2.26) 

and 

{bip),b\p)} = {d{p),d\p)} = {27Tf2u;p6{p-q). (2.27) 
We will consider moving kinks, incoming or outgoing with a velocity v, the profile of 



which are found by applying the appropriate boost to the static solution (2.11 



Mt,x,z) = tanh(7(z - zq - vt)/D), (2.28) 
and correspondingly for the fermion modes 



riit,x,z) = ^2±lTK{l[z-zo-vt])xi{l[t-v{z-zo)],x), (2.29) 

r2{t,x,z) = ^ri{t,x,z). (2.30) 
7 + 1 

2.2 The boundary and more fermion bound states 



The boundary action (2.2) leads to the boundary conditions 



dz(l)\z=o = -^J^ir -l)\z=o, (2.31) 
Hz=o = 0. (2.32) 



These were dubbed — BC boundary conditions in [11| (see also |T^). Importantly, the static 
(but not the boosted) kink profile satisfies the boundary condition, and the kink can hence 
"go into" the wall as it is reflected. Energy is conserved with the reflected kink having 
(almost) equal and opposite velocity to the incoming one, with only a small amount of 
kinetic energy lost to scalar radiation. Fig. || shows the kink profile as it hits the boundary 
(left) and the evolution of the scalar at the boundary for various choices of the incident 
velocity (right). Unsurprisingly, for larger velocity the kink penetrates further into the wall 
before bouncing back. 

When including the boundary action term ( |2.2D it is energetically favourable for the 
scalar field to take on the value 1 (rather than —1) at the boundary |11|, and we take this 
as the initial condition. This is consistent with a kink incoming from the left, which also 
has cl){z = 0) ~ 1. 

As well as the kink, the boundary also carries localised fermion modes and we perform 
the decomposition 

^pl (t, X, z) = 0, -02 {t, X, z) = Tb (t, z)x2 it,x). (2.33) 

If the kink is far away, then near the boundary we have that (^{z) = 1, and the fermion 
equation of motion reduces to 

dzJ'B{t,z)-gJ^B{t,z) = 0, (2.34) 

X2{t,x)-a^djX2{t,x) = 0. (2.35) 
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Figure 1: Left: The incoming kink profiles before, during and after a kink-boundary collision, for 
V — 0.6. Right: The central value (j){t, z = 0) for various collision velocities. 



which has the (single) solution 



(2.36) 



for the z-dependence and a left-handed chiral fermion for the boundary field. Note that 
we are forced to impose the condition ipi{z = 0) = in order to find normalisable (and 
numerically stable) solutions only. 

As a result, the boundary modes have opposite helicity to those on the kink, and in 
order to quantize the system we expand the spinor in terms of these, 



X2{t,x) 



(fp 1 



(27r)2 2u;p 

2.3 Bogoliubov coefficients 

For a static kink we have for each mode p 



e{p)e'P''U{-p) + f\p)e-'P''U{-p) 



(2.37) 



^^+{t,x,z,p) 







[cosh.{z/D)]9D I ieip^u{p) 



it,x,z,p) 



Mr 



K 



[cosh.{z / D)]9D \ie-'P''U{p) J ' 

(2.38) 



and for the boundary 

^'^+ (t , X, z, p) = Mb exp(5z) 



eip^U{-p) 



{t,x,z,p) =MBew{9z) 



-^p^U{-p) 




(2.39) 



In order to quantize the world- volume fields we need to introduce a set of creation/annihilation 
operators. We consider the system in the asymptotic past, long before the collision, and 
in the far future, long after the collision. In each of these eras the world volume fields may 
be expanded in terms of perturbations around the vacuum, as in standard quantum field 
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theory. As such we write the field operator in the distant past as 

cfp 1 



^{t,x, z) 



(27r)3 2u;p 
d?p 1 



h{p) ^>f+ {t,x,z,p) + d^ (p) (,t,x,z,p) 



(27r)3 2LUp 
+bulk continuum modes 



e{p)^Z+ it, x,z,p)+f'' {p)^f- it, x,z,p) 



(2.40) 



suitably boosted to velocity v. In this way we find that: creates positively charged, 
positive helicity particles on the kink; creates negatively charged, negative helicity anti- 
particles on the kink; creates positively charged, positive helicity particles on the bound- 
ary; p creates negatively charged, negative helicity anti-particles on the boundary. We 
shall not be concerned with the bulk contimiTim modes. 

We then let the kink collide with the boundary and move away again. The resulting 
fermion field operator can then be re-expanded (in the asymptotic future) in terms of the 
same modes, boosted in the opposite direction^ 

d'^p 



*(i,x, z) 



+ 



1 

(27r)3 2^ 
d^p 



Hp)^^tit,^, z,p) + d^p)^^Jit,x, z,p) 



1 



~eijp)^^tit,x,z,p) + f{p)^>^tit,x,z,p) 



(27r)3 2iOp 
-l-bulk continuum modes 



(2.41) 

(2.42) 
(2.43) 



By comparing the e*^'^ and U dependence of the initial and final states we find that 
the two are related through Bogoliubov coefficients a,. Pi as 



^fn^it,x,z,p) = aK+ip)K^it,x,z,p) + l3K+ip)^outit^^,z,-p) + b.c.m. 



B- 



^f^ it,x,z,p) = aK-ip)Kut it,^, z,p) + pK-ipJ'^outit^x, z, -p) + b.c.m, 



B+i 



^f+it,x,z,p) = aB+ip)'i'o^tit,x,z,p) + PB+ip)Kut it,x,z,-p)+b.c.m 



K-i 



^fn it,x,z,p) = aB-ip)'i>outit^^,z,p) + PB-ip)'i'outi't,x,z,-p) + b.c.m.. 



K+i 



(2.44) 
(2.45) 
(2.46) 
(2.47) 



where we have neglected contributions from the bulk continuum modes and any additional 
localised excited modes. Given that the mode functions are orthonormal we see that the 
coefficients are then simply the overlaps of the mode functions, for instance 

aKAp,t) = J dzd^x {^>^it,x,z,p)y ^>f^^it,x,z,p). (2.48) 



We have that 



bip) = aK+ip)bip) + PB-i-p)f\-p), 
dip) = a*K_ip)dip) + p*B^i-p)eH-p), 
e{p) = aB+ip)eip) + l3K-i-V)dH-p), 

hp) = aB-ip)fip) + Px+i-pPi-p)- 



(2.49) 
(2.50) 
(2.51) 
(2.52) 



= '^{vin), '^out ~ '^{vout), where in practice we determine Vout from the evolution of the kink. We 
find that Vout = —Vin to within 10~^. 
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The mode functions are all normalised which means that if the bulk modes are negli- 
gible Bogoliubov coefficients relating the bound-state modes form a unitary matrix. This 
implies that 

|aj^_P ~ |aB+|^, las-P - \aK+?', (2.53) 
\PkA^^\M\ \(5bA^ ^\fiK+\^- (2.54) 

The final particle number in for instance the K+ mode can then be found in terms of 
the particle numbers in the initial state 

{2TTf2u:pn'<+{p)5\p - q) = {h{p)%{q)) (2.55) 

= {(^aK+(mp) + [3B-[-p)f\-p)f {aK+{q)h{q) + PB^{-q)f\-q))) 

The expectation number of the particle operator then depends on the choice of initial 
quantum state. If we suppose that the state before the collision is in the vacuum of the 6, 
d, e, / operators then we see that the spectrum of particles after the collision is given by 

n^+(p) = \fiB-{p)\\ (2.56) 

n'^^ip) = \Pb+{p)\\ (2.57) 

n^+(p) = \f3K-{p)\\ (2.58) 

n^-(p) = \Pk+{p)?. (2.59) 

All that remains now is the calculation of the Bogoliubov coefficients, for this we need to 
solve the Dirac equation. 

2.4 Individual modes 

As the equations of motion are linear, we can consider each mode separately. For a given 
p we have 

iji{t,x,z,p) = e'P^Mt,z,p)U{p), Mt,x,z,p) = e'P^Mt,z,p)Uip). (2.60) 
Then the equations of motion give us 

^l(t, 2:,p) +z|p|Vii(t,Z,p) + (92^2(^,2:, 1?) -#(t,^)V'2(i,^;,p) =0, (2.61) 

^^2(i,2,p) - %|V'2(t,2:,15) +<92^i(t,z,^)) +5(/)(t,z)v5i(t,z,p) = 0, (2.62) 
where we have used (|2.19| ) 

p.gUip) = \p\ p.gU{p) = \p\U{p). (2.63) 

We solve the equations of motion numerically, starting from the initial condition of an 
incoming boosted kink and the appropriate fermion bound state 



= V ^ ^M^^^-i^ - zo)] ^ (2-64) 

^2 = ^^1' (2-65) 
7 + 1 
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for the kink mode and 



Vii = 0, (2.66) 
i^2=NB&M9z), (2.67) 

for the boundary mode. Then after the simulation has run its course we reexpand on the 
same modes, but boosted with the velocity of the outgoing kink as described above. This 
gives us the Bogoliubov coefficients. 

This covers the positive energy modes ( p.38D ( pl39| ). The negative energy modes 
( |2.38| )( |2l39 ) lead to eqs. ( |2.61| ), ( p. 62 ) with \p\ — which amounts to complex 



conjugating the equation of motion, so we should find that |ax+| = \I3k+\ = I /3k- I, 

\(y.B+\ = Iob-I, = We checked numerically that this is indeed the case, and 

so in the following we shall only quote the absolute values |aj^p, jasp, l/^i^P, 

3. The fermion spectrum 




Figure 2: The p-dependence of the Bogoliubov coefficients; of a mode on itself Ictif^sp, and 
transfered from one mode to another \j3K 

= 0.6, .9 = 2. 

We show in Fig. ^ the p-dependence of the various Bogoliubov coefficients. In this 
example, the incident velocity v is 0.6 and the coupling g is 2. Throughout the simulations, 
we use A = 2, Z? = I. 

If there were no bulk modes, nor any extra excited world-volume modes, we would 
have that jaA'P + 1/5^-1^ = 1 = |a_Bp + |/3Bp- Therefore, what Fig. § is telling us is that 
for pD<2 there is very little excitation of the bulk modes, but for higher momenta the 
bulk modes become more important. This is to be expected given that the mass gap for 
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the bulk fermions is g (where we have the |i;^|fettifc = 1), which for the simulation of Fig. 
|2| was g = 2. Now recall from ( 2.56 )-( ^?59|) that the /? coefficients give us the particle 
spectrum for the world-volume fields, if they are in the vacuum state before the collision. 
We see that the spectrum is peaked around zero world-volume momentum, with a width 
in momentum-space given by the inverse of the wall width, D. 

3.1 Large p 
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Figure 3: The momentum p* for which the fermions start to fall off the kink, compared to simple 
estimates, g-dependence (left) and v dependence (right). 



We can understand the behaviour in Fig. |2|by rewriting the equations of motion (2.61), 



(2.62) as second order differential equations, to find 

^i(t, z) = dl^i{t, z) - {p^ + g^cPit, z))Mt, z)+g (4>Mt, z) + d,Hi{t, z) ) , (3.1) 

4(t, Z) = dlij2{t, Z) - (/ + g^^it, z))Mt, z)-g {Hl{t, Z) + a,</)^2(t, z) 



The localisation of the fermion modes is due to the decrease in (f){z) and hence mass within 
the kink. Clearly, when p^ dominates g'^4>'^ as well as the other terms on the right hand 
side, the localisation will not be as manifest. 

We can estimate when this will happen by comparing p^ to the size of the other terms. 
In the kink background, we can use (|2.28| )( ^39| )( 2.30| ) to write 



7 + 1 



Hence, we expect the p term in (|3.1|) to dominate when 



P> 9, P> V9' 



C7 



P > V91- 



(3.2) 



(3.3) 



Fig. ^ (left) shows the p value at which the modes start to "fall off" (p*) plotted against g, 
for V = 0.6. This quantity p* is found from the simulations by observing where {axl starts 
to decrease. Keeping in mind to allow for an overall factor, the functional dependence 
of p* is best reproduced by the -1/57 line, orginating from the spatial derivative term in 
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( |3.l| ). The p* derived from the K and B modes are also roughly compatible. In Fig. y| 
(right) we fix g and show the dependence on v. Again, the behaviour of the K mode p* 
is roughly described by the ^/g7-dependence. When taking p* from the B mode, however, 
we observe a somewhat different dependence. Still, we believe we understand qualitatively 
the behaviour of the spectrum at high p. Our main focus will be the transfer between B 
and K modes, whih manifests itself at small p. 

3.1.1 Velocity dependence 




Figure 4: The velocity dependence of the l/^if^j-spectrum for 5 = 2 (left) and <? = 4 (right). 
Overlaid the u-dependence of the p — mode. 

As the incident velocity is increased, the amplitude of the (j){z = 0) dip at the collision 
(Fig. 1^, right) increases as the kink is able to penetrate deeper into the wall. This has a 
significant impact on the fermion spectrum, as can be seen in Fig. |^. For g = 2, the width 
of the distribution increases as v goes from 0.2 to 0.5, beyond which also the amplitude 



starts to grow. We note that this is simply the behaviour of the zero-mode [12, 11 1. For 
g = 4 the picture becomes somewhat more complicated, in that for v < 0.4, the distribution 
grows and widens, but then shrinks back to very small values for v > 0.4. Hence simple 
kinematics does not explain the behaviour. The dependence on v is oscillatory rather than 
monotonic. 

3.1.2 Coupling dependence 

We now fix the velocity and vary the coupling g. In [^] we found that the A; = mode 
has an oscillatory dependence on g, with amplitude and phase somewhat dependent on v. 
In Fig. ^ (left) we see the spectra for v = 0.4 at various values of g. The overlay shows the 
zero-mode oscillation, and to a good approximation, the spectra are almost Gaussian with 
the zero mode setting the amplitude. The width has a separate dependence on g. 

In Fig. ^ (right) we show the result when fixing v to 0.9. Again, we overlay the zero 
mode result, which sets the amplitude. But now the spectrum is (at least for 5 > 3) a 
non-monotonic function of p, with maximum away from p = 0. Again, the width has a 
separate dependence on g, but whereas for v = 0.4 the trend is for the width to shrink with 
amplitude, for v = 0.9 the opposite is the case. Note also the overall difference in vertical 
and horizontal scale between the two figures. 
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4. Conclusion 



Our purpose was to model the particle creation event, as observed by someone on either 
brane or boundary, due to the collision between brane and boundary. We have extended 
the work of by including fermions with momentum in the braneworld directions. 

We found that a non-trivial fermion spectrum emerges, in general strongly centered 
around p = but at high incident velocities it becomes peaked away from zero momen- 
tum. The oscillatory dependence, reported in |12, 11|, of the zero momentum coefficients 
on scalar-fermion coupling g and velocity v extends to finite p. The zero-mode sets the 
amplitude of the spectrum, with the width having its own dependence on v and g. In the 
previous study we included an excited zero momentum bound state as well as the ground 
state, there we found that this mode had little impact on the ground state fermion |11]. In 
principle there will also be excited modes with non-zero braneworld momentum. While the 
previous study reveals that such modes have little impact for low momentum, they may 
play a more significant role at higher momentum 

Generalising to more complicated scalar potentials and hence different domain wall 
profiles is straightforward, although one may run into trouble when solving analytically for 
the bound states. If so these can simply be found numerically. Such a system has been 
analyzed for the matter sector of the standard model , although braneworld gauge fields 
remain problematic. With the addition of extra fields we would expect the spectrum of 
braneworld particles to thermalize, eventually reaching an equilibrium distribution. 

A natural further extension of this work would be to include multiple scalar fields 
(complex, SU{2) doublets) and allow for C and CP symmetry breaking. This of course 
would be relevant for baryogenesis models, where branes deposit asymmetric amounts of 
matter and anti-matter on the world brane. 

In 3-|-l dimensions, similar methods may also be employed specifically for electroweak 
baryogenesis, by interpreting the kink as an advancing bubble wall and calculating reflection 
and transmission coefficients for fermions hitting or being hit by the CP-violating wall. We 
would then have to consider a finite temperature state of fermions, and possible include 
gauge fields in the background in addition to the scalar domain wall. 
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